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Abstract 

I We consider a coarse graining of NN interactions in coordinate space very much in the spirit 

of the weh known l^owk approach. To this end we sample the interaction at about the minimal de 
Broglie wavelength probed by NN scattering below pion production threshold. This amounts to 

i provide a simple delta-shell potential plus the standard OPE potential above 2 fm. The possible 

I simplifications in the Nuclear many body problem are discussed. 

1 Introduction 

[ The NN interaction provides a basic building block of atomic nuclei. A milestone in the development 

■ of the field was achieved when the Nijmegen group generated a fit via a partial wave analysis (PWA) to 

■ a set of about 4000 NN scattering data with x^/dof 1 [1] after charge dependence (CD) effects were 

■ incorporated and discarding about further 1000 of 3a inconsistent data. The analysis was carried out 
. using an energy dependent potential for which nuclear structure calculations become hard to formulate. 
! Thus, energy independent high quality potentials were produced with almost identical x^/dof [21 [3]. 
I Among them, the AV18 potential is directly useful for ab initio Monte Carlo calculations up to A = 
' 10 While all these potentials differ in their form, in the last years it has been realized that if 

■ CM momenta above A = 400MeV are explicitly integrated out, the remaining effective interaction 
. has appealing features. The so-called l^owk potentials [5] exhibit an astonishing degree of universality, 

produce a rather smooth interaction and weaken the strength of the interaction so that Hartree-Fock 
[ calculations may be reliable starting points for nuclear structure calculations. In the present talk we 
adress a suitable formulation of the problem in configuration space. 

Regarless of these successes it is to date unclear what is the impact of NN uncertainties on finite nuclei 
due to our ignorance on short distances. Relevant length scales are a) The mean interparticle separation 
distance d = 1.8fm as obtained from Nuclear matter saturation density po = ^/d^ = 0.17fm~^, b) 
The Fermi momentum kp = {3/2)3 /d ~ 250MeV which gives a wavelength of about h/kp = 0.8fm, 
c) Minimal relative CM de Broglie wavelengh corresponding to the pion production threshold A = 
h/ ^/Mmttlt^ ~ 0.5fm and d) The pion Compton wavelength l/m.,^ = 1.4fm. The situation is presented 
pictorially in Fig. [1] suggesting that for the description of the ground state in light nuclei both the short 
distance core and the role of explicit pions become marginal. This was recognized long ago [6] where the 
bulk of ^He and "^He could be described with a pionless and soft-core potential which just reproduced 
the S-wave phase shifts up to -Elab < lOOMeV. Actually, we expect this feature to hold for light nuclei. 

*Presented by R.N. P. in From Quarks and Gluons to Hadrons and Nuclei 33rd Course International School of Nuclear 
Physics, Erice, Sicily, 16-24 September 2011. 
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Figure 1: Left panel: Cartoon of a nucleus, displaying the size of the nucleons as compared 
to the typical distance to nearest neighbors and the shortest wave-lengths wave functions. 
Right panel: The AV18-potential in the ^5*0 channel. Superposed are the eigen wavefunctions 
at zero energy, Elab = 100 MeV and Elab = 350 MeV. 



2 The delta-shells (DS) potential 

The basic observation was made long ago by Aviles [7] and recently rediscovered in the context of 
renormalization of chiral forces jH]. If the two-particle CM wave numbers are limited to a range Ak 
only gross information can be determined in an interval Ar, (see e.g. Fig. [1]) with ArA/c ~ 1. Thus, for 
Afc = A ~ 400MeV we have Armin ~ 0.5fm. This uncertainty suggests that for a limited energy range 
the potential only needs to be known in a limited number of points. With this in mind we consider a 
neutron-proton (np) potential as a sum of 5 functions 

i=i 

where = Mj^l2 is the reduced np mass of the system, the Aj coefficients are strength parameters and 
Ti are the concentration radii. In that case we may determine the s-wave as 

u{r) = sm{kr + ^j+i/a) , Ti < r < r^+i , (2) 

where 6i+i/2{k) is the accumulated phase shift at the mid-point rj+1/2 and 6{k) = 5jv+i(/c). Matching 
the discontinuity of log- derivatives at r = r^, we simply get 

k cot{kri + - k cot{kri + 5i_i/2) = Aril (r^) , (3) 

where U{r) = 2^V{r) = MV{r) is the reduced potential. The regular solution at the origin requires 
Si{k) = 0. If we take the limit Ar — )• we can define 6{k,ri) = 6i{k), to get 

= -yU{R) sm\kR + 6{k, R)) + 0{Ar^) . (4) 
dR k 

which is the variable phase equation [H] up to finite grid corrections and can be interpreted as the change 
in the accumulated phase when a truncated potential of the parametric form U{r)6{R — r) is steadily 
switched on as a function of the variable R. This equation and its generalization to coupled channels 
has extensively been used to treat the renormalization problem in NN scattering in Refs. [TOl [HI [12] . 
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Figure 2: Upper panel: The realistic AV18 potential (^5'o-cliannel) and its delta-shell repre- 
sentation (left) and the coarse grained potential (right). Lower panel: Accumulated phase 
shifts for both function of the LAB energy (in MeV) indicating the number of grid 

points. Both resulting phase-shifts coincide. 



The low energy expansion of the discrete variable phase equations was used in Ref. |T3] to determine 
threshold parameters in all partial waves. The relation to the well-known Nyquist theorem of sampling 
a signal with a given bandwidth has been discussed in Ref. [8]. Of course, this DS approximation to 
the potential can be most immediately used as a numerical method to solve the scattering problem, 
which would become exact for Ar — )■ 0, if we take the weights given by the potential Aj = U{ri)Ari. 
As an illustration we show in Fig. [2] the phase-shifts obtained for the ^S'o-AV18 potential for several 
values of Ar. Convergence to the phases to four significant figures is achieved for Ar = O.Olfm. 
The equidistant discretization corresponds to the trapezoidal rule and one could improve by a more 
sophisticated method, a relevant issue when the interaction is known a priori. 



3 Coarse grained local potentials 

Another, more fruitful and economical, perspective already pursued by Aviles corresponds to consider 
the weights themselves, Aj, as fitting parameters to the phase-shifts, since anyhow the potential at short 
distances is unknown and will be determined from the data. 

If we take just one delta-shell in S-wave we may determine both the point r = r^ and its corresponding 
strength, Ac from the scattering length ao and the effective range ro, defined from 

kcotSo{k) = — - + -rok'^ + ... (5) 

For instance, for the ^5*0 case one has agiso = — 23.74fm and Tqisq = 2.77fm which yields Xisq = 
—0.4626 and r^iso = 1.99fm whereas for the ^5*1 channel one gets 00,35^ = 5.42fm and ras^ = 1.75fm 
giving A35J = —0.911 and r^^s-^ = 1.53fm. The corresponding phase shifts are reproduced to about 
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Figure 3: Phase-shifts fitted to the PWA of Nijmegen in the lowest partial waves. 



Elab ^ 50MeV. One can improve on this by including more delta-shells. A good fit to the PWA is 
obtained for (Aj in fm~^ and rj in fm) 

{Xi,ri) = (-0.568, 1.56), (-0.023, 3.47), , 

(Ai,ri) = (-0.951, 1.35), (-0.052, 2.60), . (6) 

This shows that one can consider the grid points as well as the weights as fitting parameters. The result 
for 5 equidistant points with Ar = 0.7 fm is shown in Fig. |2] (right panel). Of course, the existence of 
finite experimental errors helps in decreasing the number of coarse grained grid points. 

We have carried out preliminary fits to the NN database [1] with a pion tail with an average 

= 138MeV starting at 2fm for partial waves with J < 4 and about 40 parameters with x^/dof < 1—2. 
The result for low partial waves is shown in Fig. [31 The full PWA using a CD-OPE potential tail with 
the pertinent electromagnetic corrections to the PWA database will be presented elsewhere. 

It is straightforward to look at the deuteron by analysing the ^Si Di channel for negative energy. 
The results can be seen in Table [H The deuteron wave functions as well as the corresponding charge 
form factor is displayed in Fig. |H The peaks in the wave functions correspond to the discontinuity in 
the derivatives at the chosen grid points which, as wee can see, does not become dramatic for the form 
factor at the considered g's. 

Fourier transforming the DS potential in the S-waves gives 

{k'\V^i'\k) = Y.{\.)i,',r^3v{k'n)3i{kn) , (7) 

i 

which is a finite rank separable potential, a representation which proved very handy in the past for 
few-body and nuclear matter calculations (see e.g. [H]). 
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Figure 4: Left panel: The deuteron wave functions Right paneh Charge form factor. 



Table 1: Deuteron parameters 





Delta Shell 


Empirical 


Nijm I i 


Nijm II [2J 


Reid93 [2] 


AV18 [3] 


7(fm-^) 


0.230348 


0.231605 


Input 


Input 


Input 


Input 


V 


0.02488 


0.0256(5) 


0.02534 


0.02521 


0.02514 


0.0250 




0.8768 


0.8781(44) 


0.8841 


0.8845 


0.8853 


0.8850 


rm(fm) 


1.9676 


1.953(3) 


1.9666 


1.9675 


1.9686 


1.967 


^^(fm^) 


0.2693 


0.2859(3) 


0.2719 


0.2707 


0.2703 


0.270 


Pd 


5.498 


5.67(4) 


5.664 


5.635 


5.699 


5.76 



We show in Fig.[5]a comparison of the Fourier transformed DS potentials in the ^5*0 and '^5'i channels 
with with parameters as in Eq. to the corresponding diagonal elements of the Vio„k potentials [5], 
obtained from the AV18 [3] interaction. While the resemblance is indeed rather close, we do not expect 
a perfect description since the way the scattering problem is treated in the DS case is different as in 
the Viowk case. We have checked that one can represent quite accurately the current diagonal pieces of 
the Viowk potentials by Eq. ([7]), but does not necessarily reproduce the off-diagonal matrix elements 
constructed in the Vfowk approach from the truncated the half-off shell Lippmann-Schwinger equation 
by a specific block-transformation method. 
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Figure 5: Comparison of the coordinate space coarse-grained potentials in the ^5*0 and ^5*1 
channels with the corresponding l^iowk potentials [5J, obtained from the AV18 |P] interaction. 
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Figure 6: Binding energy of closed-shell nuclei ^He, ^^O and ^°Ca using a HO shell model 
Slater determinant as a function of the msr when the relevant phase shifts are reproduced up 
to -Elab = lOOMeV and -Elab = 350MeV. We compare to other methods (see main text). 



4 Closed-shell nuclei 

When switching from the NN problem to the many body nuclear problem the features and the form of the 
interaction are relevant in terms of computational cost and feasibility. We coarse-grain the interaction, 
but keep the exact kinetic energy, so that for two nucleons at a relative distance \xi — X2I 7^ the 
interaction vanishes, and hence the wave function becomes a Slater determinant of single particle states 

1p{Xi, . . . , Xa) = A [0ni,Zi,s,msi,t,mtj {^i) . . . 4>nA,lA,s,msA,t,mn {^a}] ■ (8) 

We use Harmonic Oscillator single particle wave functions with oscillator parameter b, where the 
spurious CM motion is exactly subtracted, for the shell-configurations ^He:(ls)^, ^^0:(ls)^(lp)^^ and 
^°Ca:(ls)^(lp)^^(2s)''(l(i)^°. Generally, for double-closed shell nuclei one has 

{V2)A=Y.9niJs{nl\V^\nl), (9) 

nlJS 

in terms of the relative matrix elements and gnisj depends on the Talmi-Moshinsky brackets 0. Using 
the single delta function which is just fixed with the S- waves scattering lengths and effective ranges (see 
below Eq. one obtains at the minimum i?(^He) = 20MeV. In common with other soft potentials [B] 
the interaction does not require strong correlations in the many-body wave function. This is due to the 
fact that since the phase-shift is reproduced to about -Elab = 50MeV the core may be ignored. Clearly, 
if we insist on reproducing up to -Elab = 350MeV a strongly repulsive DS contribution emerges and thus 
a product wave function is not appropriate. One can improve on this by adding more deltas as in Eq. 
but keeping the fit to £'lab ^ 100 MeV, in which case B{^B.e) = 24MeV. This is surprisingly close to 
the Green Function Monte Carlo (GFMC) AV18 [1] and the UCOM method [15] without three-body 
forces and complies to a cancellation between the core in the and the correlations the wave function. 

The results for the binding energy as a function of the corresponding msr radius are presented in 
Fig. El We compare with the UCOM method jl5j, Brueckner-Hartree-Fock (BHF) [16] and Coupled 
Cluster (CC) [T7]. In the UCOM method [15j a unitary local transformation generates a smooth nonlocal 
interaction from the AV18-potential while the wave functions are the same. As advertised, our results 
depend on the fitted energy range, somewhat resembling analogous ambiguities as those of the UCOM. 



^For instance, for ''He one has a m.s.r. = 35/2 V2 and for Ru{r) = 2e'fi^/^b^/ 



For an Android implementation of these calculations see e.g. http://www.ugr.es/~amaxo/eiiidroid/| 
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5 Conclusions 



We have shown how samphng of the NN interaction by a delta shell potential with a resolution de- 
termined by the deBroglie wavelength of the most energetic particle provides a coarse graining in 
configuration space, analogous to the Viowk approach. However, rather than transforming a high quality 
potential we suggest to determine the NN coarse grained interaction directly from the scattering data. 
A preliminary fit to the np phase shifts in the Nijmegen data base to all partial waves with j < 4 
requires about 40 fitting parameters yielding x^/d.o.f. < 2 (less than 1 in some waves). Deuteron 
properties show good agreement with empirical values and other calculations. Harmonic oscillator shell 
model variational calculations of nuclear binding energies provide results at the 20-30% accuracy. 

The work is supported by Spanish DGI and FEDER funds (grant FIS2008- 01143/FIS) and Junta 
de Andalucfa (grant FQM225). R.N.P. is supported by a Mexican CONACYT grant. 
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